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Abstract The large majority of EGRET point sources 
remain to this day without an identified low-energy coun- 
terpart. Whatever the nature of the EGRET unidenti- 
fied sources, faint unresolved objects of the same class 
must have a contribution to the diffuse gamma-ray back- 
ground: if most unidentified objects are extragalactic, 
faint unresolved sources of the same class contribute to 
the background, as a distinct extragalactic population; 
on the other hand, if most unidentified sources are Galac- 
tic, their counterparts in external galaxies will contribute 
to the unresolved emission from these systems. Under- 
standing this component of the gamma-ray background, 
along with other guaranteed contributions from known 
sources, is essential in any attempt to use gamma-ray ob- 
servations to constrain exotic high-energy physics. Here, 
we follow an empirical approach to estimate whether a 
potential contribution of unidentified sources to the ex- 
tragalactic gamma-ray background is likely to be impor- 
tant, and we find that it is. Additionally, we comment 
on how the anticipated GLAST measurement of the dif- 
fuse gamma-ray background will change, depending on 
the nature of the majority of these sources. 
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1 Introduction 

The EGRET telescope aboard the Compton Gamma- 
Ray Observatory detected, during its nine years of op- 
eration, not only 271 point sources [T] but also diffuse 
emission [2j|3]. Based on the origin of the diffuse pho- 
tons, we can classify this emission as Galactic if it is 
produced within the Milky Way, or as extragalactic if it 
originates from larger, cosmological distances. The two 
types of emission can be identified, to a certain extent, 
from their spatial distribution: Galactic emission is ex- 
pected to be enhanced near the Galactic plane, while 
extragalactic emission is expected to be largely isotropic 
- however, the Galactic emission does not quickly fall 
off to zero as we move away from the plane. Disentan- 
gling the Galactic and extragalactic components requires 
modeling the Galactic emission through cosmic ray prop- 
agation models (e.g. [2l[6]) and subtracting it from the 
all-sky diffuse signal. The uncertainties involved in this 
process are nontrivial. Residuals from inadequate or im- 
perfect treatment of the dominant Galactic component 
can severely contaminate the determination of the extra- 
galactic emission [HISIIB]. 

Diffuse emission can also be characterized, depending 
on the process by which it is produced, as truly diffuse 
or as unresolved point source emission. Truly diffuse is, 
for example, the emission resulting from cosmic ray in- 
teraction with interstellar matter in the Milky Wajo, as 
well as the emission from particles accelerated in shocks 
at the outskirts of cosmological large-scale structures. In 
contrast, unresolved point source emission is the emis- 
sion produced by a collection of faint, unresolved point 
sources such as blazars, perceived as diffuse due to limi- 
tations in telescope sensitivity. 

Although many sources have been suggested to be the 
origin of the diffuse emission, there are some guaranteed 

^ Note however that enhanced emission close to cosmic-ray 
sources could in principle be resolved as point source emission 
and would thus count as "unresolved point source emission" 
in this classification scheme. 
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contributions. Any known class of gamma-ray sources 
with some already identified members must have at least 
some contribution to diffuse emission. This would origi- 
nate from the collective emission from fainter members, 
unresolved by EGRET. Prime examples of such classes 
are blazars, which may be a dominant component of the 
extragalactic diffuse emission, e.g. [7l[8ll9lfT0] . as well as 
normal galaxies, e.g. [TTlfT^ . or pulsars (for the case of 
the Galactic diffuse emission), e.g. [l3 lfT4l[T5] . 

It is also possible that unresolved sources of the same 
class as unidentified EGRET sources have some appre- 
ciable contribution to the extragalactic background. Al- 
though their nature remains unknown, it is reasonable 
to believe that there is a large number of fainter, un- 
resolved objects of the same class, making some con- 
tribution to the diffuse emission. In addition, unidenti- 
fied sources are the most numerous group of gamma-ray 
sources. If they represent yet unidentified members of 
some known class of gamma-ray emitters (e.g. blazars), 
then excluding them from any calculation of the con- 
tribution of the parent class to the diffuse background 
would lead to a significantly underestimated result, due 
to an incorrect normalization of the bright-end of the 
gamma-ray luminosity function. If they represent an un- 
known class of gamma-ray emitters (which is likely - see 
e.g. [16j). then the contribution of their unresolved coun- 
terparts to the diffuse emission would significantly limit 
the diffuse flux left to be attributed to known classes and 
to truly diffuse emission. 

Similarly, investigating the diffuse emission from un- 
resolved unidentified objects is important regardless of 
the location of these objects (whether they are Galac- 
tic or extragalactic). If they are extragalactic, then un- 
resolved objects of the same class contribute to the dif- 
fuse extragalactic gamma-ray background. Alternatively, 
if they are all Galactic, objects of the same class in 
other unresolved galaxies enhance the contribution of 
their hosts to the gamma-ray backgrouncJl. 

Hence, some contribution of unresolved unidentified 
sources to the extragalactic diffuse background is certain. 
It is therefore clear that until we either (a) resolve the 
issue of the nature of unidentified sources or (b) derive 
some strong constraint indicating that a possible contri- 
bution of such unresolved objects would indeed be minor, 
we cannot hope to be confident in our understanding of 
the origin of the extragalactic gamma-ray background. 
However, predictions for the level of their collective con- 
tribution involve important uncertainties: due to lack of 
identification of low-energy counterparts, we have no es- 
timates of distance, and therefore no estimates of the 
gamma-ray luminosities of these sources. For this rea- 
son, very few constraints can be placed on their cosmic 
distribution and evolution. 



^ Additionally, unresolved sources of the same class in our 
own galaxy contribute to the diffuse emission of the Milky 
Way. 



In this work, we approach the problem from a purely 
empirical point of view. Instead of attempting to predict 
the level of a diffuse component due to unresolved objects 
of the same class as unidentified EGRET sources, we try 
to assess whether there are any quantitative indications 
that this component is, in fact, minor. Under the as- 
sumption that the majority of the unidentified EGRET 
sources can be treated as members of a single class of 
gamma-ray emitters, and for the case that this class 
consists of extragalactic objects, we try to answer the 
following two questions: (I) is it likely that unresolved 
objects of the same class could have a significant con- 
tribution to the extragalactic gamma-ray background at 
least in some energy range, and (2) how would the collec- 
tive spectrum of their emission compare to the measured 
spectrum of the extragalactic gamma-ray background 
deduced from EGRET observations. The observational 
input constraining our calculations will be the number 
distribution of unidentified sources with respect to flux, 
and the spectral index of each source. We also exam- 
ine how we expect GLAST observations to change our 
knowledge of the nature of unidentified sources, based 
on the insight gained from our analysis. 

This paper is structured as follows. In ^ we sum- 
marize the formalism used to derive the extragalactic 
gamma-ray background component due to unresolved 
unidentified sources under our set of assumptions. In Sj3] 
we describe our results and their implications. Prospects 
for the GLAST era are discussed in !|4l Finally, we con- 
clude and discuss our findings in [JS] 



2 Formalism 

In this section, we describe how we can use the flux dis- 
tribution of unidentifled sources to extract information 
about a possible contribution of unresolved sources of the 
same class to the extragalactic diffuse gamma-ray back- 
ground. The assumptions we will make to proceed are 
that (a) the unidentifled sources can indeed be viewed 
as a class of objects, so that constructing a flux distri- 
bution is meaning fuO and (b) close to the EGRET flux 
limit, the flux distribution does not evolve drastically so 
that an extrapolation of the measured flux distribution 
to lower fluxes is representative of its behavior in the 
low-flux regime. This assumption is less likely to hold as 
the limiting flux to which we are extrapolating becomes 
lower: the flux distribution will eventually exhibit a break 
due to cosmological effects and/or luminosity evolution. 

The question we will seek to answer is how far in the 
low-flux regime our extrapolation must continue before 
we get a significant contribution of unresolved unidenti- 
fied sources to the gamma-ray background. If the answer 

^ In reality, even if most of the unidentified sources do be- 
long to the same class, there is always going to be some con- 
tamination of different class objects; however, it is conceivable 
that such contamination is small. 
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Fig. 1 Flux distribution of unidentified sources. The error 
bars represent la Poisson errors. Dashed line: power-law fit 
in the range 13 x 10"* cm"^ s"^ < F < 55 x 10"* cm"^ s"\ 
Solid line: parabolic fit in the range 3 x 10"* cm~^ s"^ ^ F < 
13 X 10"*cm"^s"\ Thick solid line: Fn,in- 

is "not very far" , then we might expect that the actual 
flux distribution of the unresolved sources does indeed 
resemble our assumed form and that unresolved, uniden- 
tified sources make up a considerable fraction of the ex- 
tragalactic diffuse background. On the other hand, if we 
need to extrapolate the flux distribution down to fluxes 
very low compared to the resolved flux range, then it is 
quite unlikely that our extrapolation is valid through- 
out the flux regime we are using it. In such a case it is 
doubtful that the actual flux distribution of the unre- 
solved unidentified sources is indeed such that unidenti- 
fied sources make up a significant portion of the extra- 
galactic diffuse background. 



2.1 The cumulative flux distribution 

In our calculation, we will only include in the "uniden- 
tified object" class those 3rd EGRET catalog sources 
which still remain without a suggested low-energy coun- 
terpart candidatt0. If any sources originally included in 
the 3rd EGRET catalog listing of unidentified sources 
have been omitted when they should in fact have been 
included in our object sample, the error that this omis- 
sion incurs is consistently toward the side of underes- 
timating the importance of the contribution of uniden- 

* Sample compiled and actively maintained by C. Brown 

and available online at 

http:/ /home. uchicago.edu/~carolynb/unidentified_sources', 
it IS important to note that suggested associations with 
lower-energy counterparts included in this compilation are 
simply results reported in recent publications; no effort 
was made to evaluate the significance and validity of these 
counterparts by a single uniform standard. 



tified sources to the diffuse background. The opposite 
extreme would be if we had chosen to include all sources 
marked as unidentified in the 3rd EGRET catalog. In- 
cluding more sources in the resolved sample than we 
did, would significantly enhance the bright-end of the 
flux distribution; this would also imply more faint-end, 
unresolved, unidentifled sources, which would support a 
stronger unidentified source component in diffuse radi- 
ation. In a future version of this work, we plan to as- 
sess quantitatively the uncertainties introduced by our 
choice of the "unidentified sources" sample by examining 
to what extent our conclusions may change when differ- 
ent sets of EGRET sources, satisfying different criteria, 
are used. 

The cumulative flux distribution, iV(> F) (number 
of sources with flux above F versus F) of these sources 
is plotted in Fig. □ Here, F is the mear|3 (P1234) pho- 
ton flux in energies above 100 MeV quoted in the 3rd 
EGRET catalog [Ij. The dashed line is the power- law fit 
to the data in the fiux interval between 13x 10^* cm~^ s~^ 
and 55 x IQ-** cm^^ g-i, Nbi> F) = CF^"" , with Fg = 
F/(10-Scm-2s-i), InC = 8.32 ± 0.36 and k = 1.67 ± 
0.11. In the "faint source" end of the cumulative flux 
distribution, the data necessarily deviate from the power 
law fit due to the finite sensitivity of the telescope. How- 
ever, because of differences in sky coverage, the angu- 
lar dependence of the diffuse background, and source 
variability, there is no sharp cutoff in the fluxes of re- 
solved objects. Instead, the cumulative flux distribution 
gradually flattens. The thin solid line shows a polyno- 
mial fit to the data in the faint source end (flux inter- 
val between 3 x 10^^ cm^^ s"^ and 13 x 10"^ cm^^ 
Nf{> F) = Ao - AiFs - A2-FI , with Ao = 119.5, 
Ai = 1.837 and A2 = 0.24214. The subscripts h and 
/ in A^(> F) refer to the "bright" and "faint" end of the 
cumulative flux distribution respectively. 

If an extrapolation of the power law fit to low fluxes 
is representative of the number of existing sources in the 
faint source end, then the differential flux distribution 
(number of objects with fluxes between F and F + dF) 
of existing sources (both resolved and unresolved) of 
the "unidentified class" in units of (10~^ cm~^ s^^)^^ is 
\dN|dF\^, = KCFg^"-^. In the faint end of source fluxes, 
the number of unresolved sources with fluxes between F 
and F + dF is the number of existing sources minus the 
number of resolved sources, or 

^ dF^ dFs [nCF^--^ - - 2A2FS] . (1) 
2.2 The spectral index distribution 

The second important observational input in our calcu- 
lation is the distribution of spectral indices of unidenti- 
fied objects. In our analysis we will adopt the assumption 

^ Note that we have not tried to address strong variability 
of the sources. 
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that the spectral index distribution of unresolved uniden- 
tified sources is the same as that of the resolved uniden- 
tified sources. The latter can be deduced from measure- 
ments of the spectral index a for each of the resolved 
unidentified sources. 

Figure [2] shows a histogram of the spectral indices of 
the resolved unidentified sources (solid line). Note that 
the typical measurement uncertainty for any single spec- 
tral index (thick solid line) is comparable with the spread 
of the distribution, so that the spread of a simple binning 
of spectral indices might not in fact give us information 
about the underlying distribution of the spectral indices 
of the sources, but rather be representative of the un- 
certainty of each single measurement. This problem is 
not unique to unidentified sources. It is also a problem 
in measuring the spectral index distribution of blazars, 
where the usually derived concavity [7] of the collective 
unresolved blazar spectrum may simply be the result of 
overestimating the spread in the spectral index distribu- 
tion [13. 

Following an analysis similar to that of [17] for the 
case of blazars, we assume that the intrinsic spectral in- 
dex distribution of unidentified sources can be approxi- 
mated by a gaussian. We then use a maximum-likelihood 
analysis which takes into account the individual errors of 
measurement of a for each sourccQ to estimate the pa- 
rameters of the distribution, obtaining a mean of qq = 
2.38 ±0.03 and a standard deviation of CTq. — 0.19 ±0.03. 
The errors quoted are Icr uncertainties of each parame- 
ter with the other fixed at its maximum-likelihood value. 
The maximum-likelihood distribution is plotted with the 
dashed line in Fig. [2l and is narrower and displaced to 
lower spectral indices with respect to the histogram (a 
result reflecting the fact that sources with higher val- 
ues of a also tend to have a larger error of measurement 
in Of). However, due to the large systematics associated 
with the observational determinations of the extragalac- 
tic gamma-ray background spectrum at high energies, 
our conclusions are not very sensitive to the small dis- 
placement of the peak. 
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Fig. 2 Spectral index distribution of unidentified sources. 
Solid line: histogram of the EGRET data. Dashed line: 
maximum-likelihood gaussian. Thick solid line: typical un- 
certainty of individual spectral index determination using 
EGRET data. 

measurement of -F. The l/47r normalization factor comes 
from assuming an isotropic distribution of sources, the 
collective emission of which is uniformly distributed over 
the celestial sphere. Using Eq. ([1]) for the distribution of 
fiuxes, the maximum-likelihood gaussian of |j2T2]for the 
spectral index distribution, and integrating over the de- 
sired flux interval and over all possible spectral indices, 
we can then calculate the collective diffuse emission due 
to unresolved sources of the "unidentified class" with 
fiuxes between Fmin and -Fmax- We will take Fmax = 13 x 
10~^cm~^s~^ (the lower end of the "bright source" fiux 
interval, where we assume that all sources with fluxes 
above this limit have been resolved by EGRET), i^min 
(the low-flux limit of the flux-function extrapolation) is 
the free parameter in our calculation. 



3 Results 



2.3 Contribution to the extragalactic gamma-ray 
background 

The extragalactic gamma-ray background is described 
by the differential photon intensity Ie (photons per unit 
area-time-energy-solid angle) . Each unresolved source of 
flux F has a contribution Ie,i{F) to the diffuse emission 
which is given by 

where a is the spectral index of the source, and Eq — 
100 MeV is the lowest photon energy included in the 

^ We do this by introducing the true spectral indices of the 
sources as nuisance parameters and by marginalizing over 
them. 



Integrating Eq. ^ we can now calculate the contribution 
of unresolved unidentified sources to the diffuse back- 
ground as a function of Fmin. Note that the overall level 
of the associated diffuse emission only depends on where 
the extrapolated power law breaks. The thick solid line 
in Fig. [1] shows the value of this cutoff so that the con- 
tribution of unidentified sources does not overtake the 
original EGRET estimate of the extragalactic gamma- 
ray background [3] in the range 100-300 MeV plus sta- 
tistical erroi0. We would only need extrapolate the cu- 
mulative fiux distribution for slightly more than an order 

^ This is a conservative choice, since all more recent esti- 
mates of the extragalactic diffuse background give an EGRB 
level lower than the original estimate of '3". If we had instead 
adopted a lower EGRB intensity, the thick solid line of Fig. [1] 
would move to higher fluxes, and this would make an impor- 
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Fig. 3 Dashed line: Sreekumar et al. (1998) determination of 
the EGRB. Dotted hue: Strong et al. (2004) determination of 
the EGRB (best guess). Solid lines: Strong et al. (2004) deter- 
mination of the EGRB (systematics-based limits) . Thick solid 
line: collective spectrum of unresolved unidentified sources 
(this work). 



of magnitude below the lower limit of the resolved flux 
range to have unresolved unidentified sources comprise 
most of the extragalactic diffuse background, at least at 
low energies. This is not an extreme extrapolation, and 
therefore a significant contribution by the "unidentified" 
class to the diffuse background is likely. 

Figure [3] shows the cumulative emission spectrum 
of unresolved unidentified sources, overplotted with the 
spectrum of the extragalactic diffuse emission derived 
from EGRET observations. The dashed line shows a single- 
power-law fit to the Sreekumar et al. (1998) determina- 
tion of the extragalactic gamma-ray background (EGRB) 
[3] . The dotted line is the more recent redetermination of 
the gamma-ray background by Strong et al (2004) |6J, in 
which they used their more detailed model of the Milky 
Way diffuse emission to subtract the Galactic component 
from the EGRET diffuse sky map. The solid lines are the 
systematic uncertainties in the EGRB determination of 
[6] , entering through their model of the Galaxy. Our cal- 
culation of the spectrum of the unresolved unidentified 
source component is shown with the thick solid line. At 
low energies, where the systematics are low, the uniden- 
tified component spectrum is in excellent agreement with 
the EGRB observational spectrum of At higher en- 
ergies, where the systematics are large, the unidentified 
component spectrum is largely within systematics except 
at very high energies. If unidentified sources are indeed 
a dominant contribution at relatively low energies, then 
this result may be perceived as a tantalizing hint that at 
the highest energies of the EGRET range a new type of 
contribution, (e.g. from annihilating dark matter) may 
start to be important at a few tens of GeV. 

tant unidentified source component of the EGRB even more 
likely. 



4 Prospects for the GLAST era 

The launch of GLAST in 2007 wiU provide us with signif- 
icant new insight about the nature of unidentified sources 
and their possible contribution to the extragalactic dif- 
fuse background. The ideal solution to the unidentified 
source puzzle would be, of course, the direct positional 
association of all unidentified sources with undisputed 
low-energy counterparts. This would then allow us to 
build more confident models for the unresolved members 
of these classes of objects. However, such an outcome is 
unlikely, as the large number of possible counterparts and 
the large number of sources which we expect GLAST will 
be able to resolve make multi-wavelength campaigns for 
every single source impractical. 

However, there is another definitive test that GLAST 
will be able to perform, which does not require confi- 
dent identification of each source to provide information 
about the likely nature of unidentified sources as a pop- 
ulation. With the increased flux sensitivity of GLAST, 
many more objects of the same class will be resolved. 
If these objects are mostly extragalactic, there will be 
an associated decrease of the extragalactic gamma-ray 
background from its EGRET levels, equal to the all-sky- 
averaged intensity of the newly resolved objects. The flux 
sensitivity of GLAST is expected to be about 50 times 
better than that of EGREtI. Therefore GLAST wiU be 
able to probe the flux distribution of unidentified sources 
down to flirxes close to i^min and definitively test our em- 
pirical estimate. 

If, on the other hand, these objects are mostly Galac- 
tic, then there will be an associated reduction of the 
" Milky Way diffuse emission rather than of the isotropic 
background. For a discussion on a possible contribution 
of a large number of Galactic point sources to the Galac- 
tic diffuse emission and the role of such sources in ex- 
plaining at least in part the origin of the GeV excess see 



5 Discussion 

In this work, we have used a purely empirical model 
to explore the possibility that unresolved gamma-ray 
sources of the same class as unidentified EGRET sources 
have an appreciable contribution to the extragalactic 
gamma-ray background. We have argued that some uniden- 
tified source contribution to the gamma-ray background 
is guaranteed. We have additionally found that (1) if 
most unidentified sources are assumed to be extragalac- 
tic, we would only need the observed cumulative fiux 
distribution of unidentified sources to extend without 
a break for a little more that one order of magnitude 
toward lower fiuxes in order to have a very significant 
contribution to the gamma-ray background - at least at 

* htt p: / /www-glcist .Stanford . edu / mission .html 
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the lower part of the EGRET energy range; and (2) the 
spectrum of the cumulative emission of such unresolved 
sources would be very consistent with the observational 
determination of ^ of the gamma-ray background from 
EGRET data, within systematics. 

We have learned that any model of the extragalac- 
tic gamma-ray background would be incomplete without 
some treatment of the unidentified source contribution. 
The results of our empirical model therefore motivate us 
to pursue specific population models for the unidentified 
sources. Although such models involve a more restric- 
tive set of assumptions and increased uncertainties, they 
can provide more concrete predictions for the luminosity 
function of unresolved objects. Additionally, if we were to 
assume that the majority of unidentified sources are in- 
deed members of a single, extragalactic class of gamma- 
ray emitters, and that unresolved members of this class 
do indeed contribute most of the extragalactic diffuse 
emission, then we can use simple evolution models to 
place limits on the redshifts of unidentified sources. We 
will pursue such models and calculations in an upcoming 
publication. 

In this work we have tried, where possible, to make 
assumptions, which, if anything, underestimate the pos- 
sible contribution of unresolved unidentified sources to 
the extragalactic diffuse background. However, some of 
our necessary working assumptions have the potential to 
overestimate the unidentified class contribution. 

First of all, we have assumed that the majority of 
the sources belong in a single class. It is conceivable that 
instead, the resolved unidentified sources are a collec- 
tion of members of several known and unknown classes 
of gamma-ray emitters. In this case, it is still likely that 
the summed contribution of unresolved members of all 
parent classes to the diffuse background is significant. 
However, the construction of a single cumulative flux dis- 
tribution from all sources and its extrapolation to lower 
fluxes is no longer an indicative test for the importance 
of such a contribution. 

Second, we have assumed that the majority of the 
sources we have used are extragalactic. It is not at all 
certain that this is indeed the case. The unidentified 
sources from the 3rd EGRET catalog [T] exhibit a strong 
concentration along the Galactic plane, and hence a sig- 
nificant fraction of the unidentified sources in the 3rd 
EGRET catalog are most likely Galactic. In the sample 
we have used, this feature is less pronounced (mainly be- 
cause recent suggestions of possible low-energy counter- 
parts refer mostly to Galactic plane objects). As a result, 
the possibility that most of the sources we have used in 
our calculations are extragalactic cannot be excluded on 
isotropy grounds. However, a mostly Galactic population 
originating in the disk and bulge is also consistent with 
constraints from the gamma-ray emission from nearby 
galaxies 11 9j . In this case the total luminosity of resolved 
unidentified sources is already a signiflcant fraction of the 
total diffuse luminosity of normal galaxies. Consequently 



the luminosity of each individual normal galaxy is signifl- 
cantly enhanced and the unidentifled sources have a con- 
tribution to the extragalactic background through their 
hosts. Ultimately, the question of whether most unidenti- 
fied sources are Galactic or extragalactic will be decided 
by GLAST. 
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